Optical Spectra and Localization of Excitons in Inhomogeneous 
Helical Cylindrical Aggregates 



Catalin Didraga and Jasper Knoester* 
Institute for Theoretical Physics and Materials Science Centre, 
University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands 
(Dated: February 2, 2008) 

We study the linear optical properties of helical cylindrical molecular aggregates accounting for 
the effects of static diagonal disorder. Absorption, linear dichroism, and circular dichroism spectra 
are presented, calculated using brute force numerical simulations and a modified version of the 
coherent potential approximation that accounts for finite size effects. Excellent agreement between 
both approaches is found. It is also shown that the inclusion of disorder results in a better agreement 
between calculated and measured spectra for the chlorosomes of green bacteria as compared to our 
previous report, where we restricted ourselves to homogeneous cylinders [J. Phys. Chem. B 106, 
11474 (2002)]. We also investigate the localization properties of the excitons responsible for the 
optical response. By analyzing an autocorrelation function of the exciton wave function, we find a 
strongly anisotropic localization behavior, closely following the properties of chiral wave functions 
which previously have been found for homogenoeus helical cylinders [J. Chem. Phys. 121, 946 
(2004)]. It is shown that the circular dichroism spectrum may still show a strong dependence on 
the cylinder length, even when the exciton wave function is localized in a region small compared to 
the cylinder's size. 
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I. INTRODUCTION 

The optical properties and optical dynamics of molecu- 
lar aggregates with a cylindrical geometry currently draw 
considerable attention. Both natural and synthetic forms 
of such molecular nanotubes are investigated. Among the 
natural systems, the rod shaped light-harvesting com- 
plexes in the chlorosomes of green bacteria are well- 
known examples i The chlorosomes of Chlorofiexus au- 
rantiacus contain tens of thousands of bacteriochloro- 
phyll molecules self-assembled in cylindrical structures 
with a monolayer wall of roughly 5 nm diameter and 
a length of hundreds of nanometers i^i^*^ The light- 
harvesting system of the bacterium Chlorobium tepidum 
also contains cylindrical aggregates, with a bilayer wall 
and a diameter of roughly 10 nm^ These natural sys- 
tems should be referred to as J aggregates, as the ab- 
sorption spectrum is red-shifted relative to the transition 
frequency of a single bactcriochlorophyll molecule. 

J aggregates with a cylindrical geometry have 
recently also been prepared via synthetic routes. 
In particular a class of substituted 5, 5', 6,6'- 
tetrachlorobenzimidacarbocyanine dyes has been 
created that forms such aggregates; the cylindrical 
geometry was revealed using cryo-TEMA^ It has been 
demonstrated that the precise morphology as well as 
the details of the optical properties depends on the 
nature of the substituents and the solvent AiS These 
synthetic cylinders usually have bilayer walls with an 
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outer diameter of about 15 nm and a wall thickness of 
4 nm. The cylinder length extends to several hundreds of 
nanometers. Recently, the (helical) arrangement of the 
molecules inside the bilayer wall of such aggregates has 
been determined for the first time. This was done for the 
aggregates of the dye 3,3'-bis(3-sulfopropyl)-5,5',6,6'- 
tetrachloro-1 , 1 '-dioctylbenzimidacarbocyaninc (C8S3) , 
by modeling the data from cryo-TEM, absorption and 
linear dichroism measurements 4fi A Frenkel exciton 
model based on two weakly interacting bricklayer mono- 
layers wrapped on cylindrical surfaces of appropriate 
diameter gave a good fit to experiment. Very recently, 
it has been discovered that bilayer molecular nanotubes 
(diameter 14 nm) may also be formed through self- 
assembly of amphiphilic hexi-peri-hexabenzocoronene 
moleculesiii Another class of synthetic cylindrical 
aggregates are those formed through self-assembly of the 
dye meso-tetra(4-sulfonatophenyl)porphyrin (TTPS4) 
in acidic aqueous solution^Sii^ From small angle x-ray 
scattering (SAXS)i^ and atomic force microscopy 
(AFM)i^ it has been concluded that these aggregates 
are hollow monolayer tubes with a diameter of about 20 
nm. Finally, it has been shown that under the influence 
of tetrahedral chemical defects conjugated polymers may 
also adopt ordered cylindrical conformations^ 

In previous model studies of cylindrical aggregates, 
we have mainly restricted ourselves to homogeneous ag- 
gregates, i.e., we have ignored the role of disorder on 
the Frenkel exciton states that determine the optical 
response li^iiSiiLiS An important simplification occurs 
in this case, as the cylindrical symmetry may then be 
used to distinguish excitons in classes (bands) of differ- 
ent transverse quantum number ^2, which describes the 
Bloch nature of the exciton wave function in the direc- 
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tion along the circumference of the cyhnderiii The ex- 
citon eigenstates may then be determined from a set of 
one-dimensional effective Hamiltonians for each value of 
k2 ■ Importantly, the introduction of the transverse quan- 
tum number dictates simple selection rules, which state 
that only states with /s2 = or fc2 = ±1 can be observed 
in linear optics. States in the k2 = band give rise to ab- 
sorption polarized along the cylinder's axis, while states 
with k2 ~ ±1 yield a polarization perpendicular to this 
axis. For each of these three bands only a few strongly 
allowed (superradiant) exciton states occur liii^ 

In spite of the fact that usually appreciable energetic 
or interaction disorder occur in self-assembled molecular 
aggregates, it turns out that a homogeneous aggregate 
model, with the simple selection rules discussed above, 
does describe the salient features of experimentally ob- 
served spectra. For instance, using this model, the ex- 
perimentally observed variation in the CD spectra of the 
chlorosomes of Chloroflexus aurantiacus was cxplainedfii 
Similarly, the polarization dependent spectra of the bi- 
layer C8S3 cylinders are rather well described neglect- 
ing the role of disorder liS We have shown, however, that 
the inclusion of static energy disorder does improve the 
comparison to experiment for the latter casein Simi- 
larly, model results for chlorosomes that do account for 
disorder also show a better quantitative comparison to 
experiment^ 

The aim of this paper is to systematically study the 
effect of static diagonal disorder on the optical spectra of 
cylindrical molecular aggregates and to study the exciton 
localization properties caused by the presence of disor- 
der. While the formalism used applies to general cylin- 
drical aggregates, we will in explicit calculations restrict 
ourselves to cylinders of the structure of chlorosomes of 
Chloroflexus aurantiacus. As observable quantities, we 
will focus on the absorption, linear dichroism (LD), and 
circular dichroism (CD) spectra. The disorder averages 
of these specta are calculated using numerical simulations 
as well as a modified version of the coherent potential ap- 
proximation (CPA) that accounts for finite-size effects. 

Localization properties have been extensively studied 
for Frenkel excitons in linear^°i^^i^^i^'^ and circular^'^i^^i^^ 
aggregates. For the higher dimensional cylindrical ag- 
gregates, such studies have not been performed yet. We 
will fill this gap by investigating the inverse participation 
ratio and an autocorrelation function of the wave func- 
tion. Using the latter, we will show a strong anisotropy 
in the localization properties of the excitons, which may 
be traced back to the recently discovered chiral behav- 
ior of the exciton wave functions on finite homogeneous 
cylindrical aggregates^^ dictated by the behavior of the 
dipolar excitation exchange interactions in the system. 
We will also demonstrate that the CD spectrum may ex- 
hibit a dependence on the cylinder's length, even when 
the excitons arc localized on regions small compared to 
the its size. 

The outline of this paper is as follows. In Scc.|n]we will 
present the exciton model and give the general expres- 



sions for the quantities of interest, in particular for the 
spectra and the localization characteristics. Section UTTI 
is dedicated to explain the modified version of the CPA. 
Results for the spectra are presented and discussed in 
Sec. IIVI while in Sec. we do the same for the localiza- 
tion characteristics. Finally, wc present our conclusions 
in Sec. ED 

II. MODEL AND QUANTITIES OF INTEREST 

a. Aggregate structure, Hamiltonian, and eigenstates. 
For the structure of the cylindrical aggregate we will use 
the same model and notation as described in Ref. ^3 and 
depicted in Figure 1. The aggregate consists of a stack of 
A^i rings of radius R (labeled ni = 1, . . . ,Ni), each con- 
taining N2 molecules. The distance between neighboring 
rings is denoted h. Neighboring rings are rotated relative 
to each other over a helical angle 7, with < 7 < 2tt/N2. 
Connecting the closest molecules on neighboring rings, 
one observes that the aggregate may be viewed as A'2 he- 
lices each containing A^i molecules, winding around the 
cylinder (the dashed line in Figure 1 indicates one such 
helix). Each molecule may now be labeled n = (711,71,2), 
where ui denotes the ring on which the molecule re- 
sides, while 712 denotes the helix on which it lies. The 
total number of molecules in the aggregate is denoted 
N = N1N2. In Ref. we explained the general nature 
of this structural model. 

All molecules are modeled as identical two-level sys- 
tems, with transition dipolcs fi^i that are equal in mag- 
nitude (/i) and have equal orientations relative to the 
frame of the cylinder at the position of the molecule. In 
particular, all molecular dipoles make an angle f3 with 
the cylinder axis (referred to as the z axis), while the 
projection of each dipolc on the a;?/-plane makes an angle 
a with the local tangent to the rings. Explicitly, the x, 
y, and z components of the molecular position vectors 
and dipole moments are given by the three-dimensional 
vectors: 

Tn = (i?cos(7i202 + ni7), J? sin(n2(/)2 + ni^),nih) (1) 
and 

fJ^n = {-fisinf3sin{n2(j)2 + ni-f — a), 

/isin/3cos(n2(/)2 + ?^i7 ~ a), Mcos/?) , (2) 

respectively, with 02 = 2tt/N2. 

The electronically excited states of the aggregate are 
described by the Frenkel exciton Hamiltonian with static 
diagonal disorder. Setting ?7, = 1. we have 

H = J2i^o + en)bibn + J2'j{n-m)bib^, (3) 

n n,m 

where bl^ and &„ denote the Pauli operators for cre- 
ation and annihilation of an excitation on molecule n, 
rcspectivelyi2L22422i2£ Furthermore, uiq is the average 
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na = 1,2,..., A^2 

FIG. 1: Cylindrical aggregate consisting of a stack of A^i rings, 
labeled ni — l,2,...,Ni, that each contain N2 molecules, 
labeled 712 = 1,2,..., N2. The arrows indicate the transition 
dipoles, which are equal in magnitude (/x) and make an angle 
f} with the cylinder axis. The projection of each dipole on the 
plane of the rings makes an angle a with the local tangent 
to the ring (see projection of one ring displayed to the right). 
Each ring is rotated with respect to the previous one over an 
angle 7, so that we may view the aggregate as a collection of 
parallel helices on the cylinder's surface. One such helix 
is indicated by the dashed curve. The label 712 in fact labels 
the helices. 

molecular transition frequency and is the static ran- 
dom energy offset at site n, induced by slow solvent ef- 
fects. J(n — m) is the excitation transfer interaction be- 
tween molecules n and m. Due to the symmetry of the 
system the interaction only depends on the relative posi- 
tions of the two molecules. The prime on the summation 
indicates that the term with n = m is excluded from 
the summation. We assume that J(n — m) results from 
dipole-dipole interactions, giving it the explicit form 

J(n - m) = ^ - 3 ■ ^ , (4) 

with r„m = r„ - Tm- 

The Hamiltonian Eq. ^ differs from the one discussed 
in Ref . ^3 only in the inclusion of random energy offsets 
£„• We will assume that the energy offsets on different 
molecules are uncorrelated and follow a Gaussian distri- 
bution V{en) with standard deviation a. Hence, each en 
is taken independently from the distribution 

To describe the linear optical response of the aggregate, 
it suffices to consider the space of onc-exciton states, i.e., 
those states in which the molecules of the cylinder share 



one excitation. The general form of these eigcnstates 
reads 

\q)=Y,V,{n)bi\g), (6) 

n 

where \g) denotes the overall ground state, in which all 
molecules are in their ground state. We have used the la- 
bel q to distinguish the N one-exciton states; the quantity 
fq{n) denotes the amplitude of the qth state on molecule 
n. These amplitudes are obtained by diagonalizing the 
N X N one-exciton Hamiltonian, which has the quantities 
ujQ + £„ as diagonal elements and the J(n — m) as off- 
diagonal ones. The eigenvectors ipq^n) will be assumed 
to be normalized to unity. Unless stated otherwise, we 
will impose open boundary conditions along the z axis. 

In the absence of disorder (tr = 0), the label q may be 
replaced by a two-dimensional label k = (fci, fc2), where 
k2 denotes the wave number describing the Bloch mo- 
mentum of the exciton state along the ring direction, 
while ki labels the Ni possible exciton states in each 
of the N2 different bands characterized by one value of 
k2 liLi^ The diagonalization then separates into N2 inde- 
pendent effective one-dimensional problems. In the cur- 
rent general case of disorder, we cannot make this decom- 
position of the quantum labels, and we will keep the gen- 
eral label q. We note that with the breakdown of the sep- 
aration into transverse (fc2) and longitudinal (fci) quan- 
tum numbers, also the strict distinction between states 
polarized parallel (fc2 = 0) and perpendicular (k2 = ±1) 
to the z axis breaks down. In general, the dipoles of 
the exciton states in disordered cylinders may have any 
orientation relative to the z axis. 

b. Optical spectra. We will be interested in calculat- 
ing the absorption, LD, and CD spectra in the presence of 
diagonal disorder. The general expressions for these spec- 
tra in terms of the one-exciton energies Eq and eigenvec- 
tor components ipq (n) follow from linear response theory 
and take the generic formii 

5(a;)= ^^X,<5(u;-i?,)y (7) 
with strength 

X, = ^^,(n)^;(m)X„,„,. (8) 

n,m 

Here, ^(a;) stands for A{uj), LD{uj), and CD(uj) in case 
of the absorption, LD, and CD spectrum, respectively, 
and the angular brackets (...) denote the average over 
the random energy offsets. Furthermore, the quantities 
^n,m are the corresponding strengths in the site repre- 
sentation, which take the form 

Xn,m = cos^ /? + ^os ^ sin^ /3 , (9) 
^n,m = cos^ /? - i/i^ COS ^ siu^ (3, (10) 
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and 



X 



CD 



COS ^) sin (2/3) cos a 
(rii — TOi)/isin^sin^ /?] . 



(11) 



for the three spectra considered. Here, ^ = [(^2— "^2)02 + 
(?ii — TOi)7] and A denotes the wavelength of the Ught. In 
the above expressions for the absorption and CD spectra, 
an isotropic average over orientations of the cyhnder has 
been used (appropriate for an isotropic solution), while 
in case of the LD spectrum, a uniform average over an- 
gles of rotation around the cylinder's axis was inferred 
(appropriate for samples with perfect alignment of the 
cylinders). Any (unlikely) correlation between orienta- 
tion and disorder realization has been neglected. For fu- 
ture reference it is useful to note that also the density of 
states, p{i^), follows the generic expression Eq. (|7I), with 
^n,m = ^n,m, implying XP^l. 

The disorder averaged spectra and density of states in 
principle may be calculated using straightforward numer- 
ical simulations, in which one generates a number of ran- 
dom disorder realizations {en} and for each realization 
performs an explicit diagonalization of the one-exciton 
Hamiltonian to obtain the quantities Eq and </Jg(n). For 
one-dimensional (linear or circular) molecular aggregates 
this is common practice. For higher-dimensional aggre- 
gates, however, such as planes or cylinders, this method 
may be quite demanding, due to the large number of 
molecules involved in these systems and the tendency 
for the collective excitations to be more delocalized in 
higher dimensions. This motivates the use of alterna- 
tive, albeit approximate, methods to calculate the aver- 
age spectra. A well-known example is the coherent po- 
tential approximation (CPA);2i4S which previously has 
been applied with success to calculate spectra of, for in- 
stance, isotopically mixed aromatic crystals (dichotomic 
disorder) one-component systems with Gaussian di- 
agonal disorder^^S*^ and two-component systems with bi- 
Gaussian diagonal disorder 

In this paper, we will use both brute-force numerical 
simulations and the CPA to calculate the spectra and 
the density of states. We will extend the usual CPA 
to account for finite-size effects (Section and show 
that this modified CPA gives excellent agreement with 
numerical simulations, in fact significantly better than 
its conventional implementation. 

c. Localization characteristics. As is well known, 
the presence of diagonal disorder leads to the localiza- 
tion of the excitonic eigenstatesi^ In order to have quan- 
titative information on this localization, we will ana- 
lyze the inverse participation ratio and an autocorrela- 
tion function of the wave function. Both quantities can 
only be addressed within a numerical simulation. The 
energy dependent inverse participation ratio is defined 
througha2ii22*^ 



E 



/P{^) ■ (12) 



The participation ratio, , is generally accepted as a 
typical value for the number of molecules participating in 
the eigenstates at energy uj. For example, a state local- 
ized on a single molecule has C{ijj) = 1, whereas for the 
completely delocalized states on a homogeneous cylinder 
C{uj) ^ 1/N . In the latter case, the precise value depends 
on the boundary conditions and on the possible combi- 
nation of degenerate complex eigenstates to real ones. If 
we use periodic boundary conditions along the z axis and 
use real sin and cos forms for the transverse and longi- 
tudinal Bloch wave functions, all states (except a few^) 
have a participation ratio 9/(4A^). For the same bound- 
ary conditions using complex wave functions, all states 
have C{lo) = 1/A^. If we use open boundary conditions 
with real transverse wave functions and accept the ansatz 
solutions for the longitudinal wave functions analyzed in 
Rcf.0, we obtain (again with a few exceptions) the value 
9/[4(/Vi-t-l)7V2] «9/(4A^). 

Alternative measures for the extent of the exciton wave 
function have also been considered, for instance, auto- 
correlation functions of the wave function were used to 
study the ring-shaped LH2 antenna systemi^^i^ Correla- 
tion functions are particularly useful when dealing with 
anisotropic higher-dimensional systems, such as cylindri- 
cal aggregates, as they allow for a study of the localiza- 
tion properties along different spatial directions. Such 
information cannot be extracted from the inverse par- 
ticipation ratio. Hence, we define the autocorrelation 
function 



C(n;tj) = /^^|v?,(m)^;(m + n)|5(w-i;,)\/p(w), 

\ q m I 

(13) 

where the summation over m extends over mi = 
1, . . . , A^i — ?7i and TO2 = 1, . . . , iV2 in order to be consis- 
tent with open boundary conditions in the ni direction. 

The generic form of C(n;Li;) on the {ni,n2) plane is a 
structure that peaks at the origin (0, 0) (where it has the 
value unity) and which (after averaging over a sufficient 
number of disorder realizations) has inversion symme- 
try with respect to the origin. The form of the peak 
shows in what direction the exciton wave functions at 
energy uj are most localized or extended. Finally, from 
the autocorrelation one can define a localization mea- 
sure (alternative to C^^iuj)) for the total number of 
molecules participating in the typical wave function at 
energy lo. We will denote this measure by N^^i{lo) and 
define it as the total number of n values (molecules) with 
C(n; cj) > C{0\Lo)/e = 1/e (e the base of the natural loga- 
rithm). It should be stressed that both the participation 
ratio and N^^^{lo) represent typical numbers and cannot 
be expected to give exactly the same result. However, 
one does expect these two measures to vary in a similar 
way with energy, disorder strength, or system size. We 
will come back to this in Sec. 
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III. MODIFIED COHERENT POTENTIAL 
APPROXIMATION 

In this section, we address some essential technical as- 
pects to use the CPA when calculating the optical spec- 
tra. The method as such is well-documented in text 
books^2iii which is why we only focus on two aspects 
that are specific to our application. The first one con- 
cerns reducing the general expression Eq. Q for the spec- 
tra to a form that can be addressed within the CPA. 
The second one concerns the treatment of finite systems, 
where periodic boundary conditions should be avoided. 

The CPA is a method that yields an approximate form 
for the disorder averaged (retarded) Green's function 

G{u;)), with 



G{lu) ^ {uji^H + ii-iY 



(14) 



Here, 1 is the unit operator, 77 is a positive infinitesimal 
constant, and H is the Hamiltonian Eq. (Q. Using the 
Green's function, we may rewrite Eq. Q for the spectra 
as 



5(c.) = -ilm{^A,(g|G(..)|g) 



(15) 



As the eigenstates \q) and the strengths Xq depend on 
the disorder realization, some care is needed to reduce 

Eq. H15|) to a form that only contains (g{uj)'^. 

We first use Eq. © for Xq in terms of Xn,m- Realizing 
that ^q{ri) = {n\q) and using the fact that the Xn m do 
not depend on the disorder realization [see Eqs. lO, (|10|1 . 
and Hll|l ] , we may rewrite 



-ilm^X„,„(n|(G(^))|m; 



(16) 



We now change to the basis of eigenstates of the system in 
the absence of disorder (cr = 0), in which case we replace 
the state label q by the two-dimensional label k = (fci , ^2) 
(eigenvectors V'k(n) and energies E]^), as explained in 
Sec.m The advantage to do this will become clear below. 
On this basis, Eq. (|16(l takes the form 



S{uj) = --ImY^X^M'm (gH) |k' 



k,k' 



where we defined 



(17) 



(18) 



We now take advantage of the fact that within the CPA 

the averaged Green's function (^G'(tj)^ is replaced by the 

Green's function of the same system in the absence of 
disorder, but with an uj dependent and k independent 
complex self-energy S(a') added to the exciton energies 



Ek. Thus, within the CPA (^G{lo)J is by definition diag- 
onal in the k basis. Hence, 

S{uj) =-ilm^Xk,k(k|(G(t^))|k) 



TT ^ ^ 



(19) 



We proceed by using the Bloch nature of the eigenvectors 
(/5k (n) in the ring direction, 

</Jk(n) = (A^2)"'/'exp[z2^fc2n2/iV2](/^fci(ni;fc2) (20) 

with the integer ^2 the transverse wave number and 
(711; ^2) the longitudinal wave functioniii Then the 
strengths ^kk' "'^kk' ^'^'^ "^k^k ^o"^ three spectra 
considered can be expressed in terms of the^tj(ni; ^2). 
This algebra was already performed in Ref. llTl yielding 
the oscillator strength Ok, the LD strength Lk, and the 
rotational strength i?k- Combining all these expressions, 
we recover for the final CPA results Eqs. (20), (24), and 
(31) of Ref. for the absorption, the LD, and the CD 
spectrum, respectively, except that 5(uj — E'k) in these 
expressions is replaced by — ilm(cj — Ey — 2(0;) -I- ir])~^ . 
In particular we find that, like in the homogeneous case, 
only terms with fc2 = 0,±1 contribute to the three spec- 
tra, with the fc2 = states having a transition dipole 
along the z axis and the other two (degenerate) bands 
having a dipole perpendicular to it. In contrast to the 
homogeneous case, however, to calculate the CPA spec- 
tra we still need the eigenenergies of all states in all ^2 
bands. The reason is that these energies occur in the 
quantity go(^) (Eq. (|23|l '). which is needed to solve for 
the self-energy. 

We finally mention that within the CPA, the density 
of states is simply given by Eq. H19|l with X^ ^ replaced 
byl. 

Since the numerical diagonalization of the N2 effec- 
tive one-dimensional problems that yield the longitu- 
dinal eigenfunctions (pkiin-i]k2) and the energies £^k is 
straightforward, the only remaining problem is to deter- 
mine the self-energy E(w). Within the CPA, S(w) obeys 
a self-consistency equation, derived by neglecting corre- 
lations between scattering events on different molecules. 
The resulting self-consistency equation reads^S*^ 



e„ - S](w) 



l-(en-E(c^))(G(c^)), 



= 0, (21) 



where the disorder average now is a simple integration 
over en, weighted by the distribution V^luo). The site 
diagonal element of the averaged Green's function is 



G(^)) =^2- 



l¥'k(n)| 



nn UJ — E'k — S(cl') + irj 



(22) 



The CPA has been developed for large systems, where 
periodic boundary conditions can safely be applied to 
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obtain the homogeneous solutions. In that case we have 
|(^k(n)|^ = with N the total number of molecules 

in the system. Substituting this in Eq. (|^ . we ob- 
serve that the diagonal element of the Green's function 
becomes site-independent, which in fact is necessary in 
order for the self-energy following from Eq. (|21|l to be 
site- independent. The latter fact is implicit in the CPA 
by assuming from the very beginning that the self-energy 
does not depend on k. Conversely, if in Eq. (|22|) we would 
boldly substitute the solution for the homogeneous solu- 
tion with open boundary conditions applied along the 
ni direction, we would obtain a site-dependent (ni de- 
pendent) self-energy, which would be inconsistent in the 
context of the CPA. 

Yet, as we will be interested in studying the effect of 
the cylinder length on the spectra, wc prefer not to use 
periodic boundary conditions in the rii direction. One 
way out of this impasse is to construct artificially a diag- 
onal element of the Green's function which is site invari- 
ant. This can be done by approximating (G(ti^))nn by its 



n 

4 E ? w N , ■ = .9oH, 



which implies in particular that the local density of states 
is approximated by the normalized total one. We note 
that in spite of this approximation, the self-energy still 
contains information about the system's finite size, as 
the finite-size energies E]^ will be used when evaluating 
go{io). Moreover, when calculating the spectra and den- 
sity of states (Eq. H19|l with the proper Xk,k) again we 
will use the exciton energies i?k as well as the strength 
Xk,k calculated for the finite homogeneous system. Thus, 
one may hope that the approximation Eq. (|23|l only af- 
fects the results for the spectra and the density of states 
in a rather weak way. In Sec. II VI we will check the validity 
of the approximation by comparing this new application 
of the CPA directly to exact numerical simulations as 
well as to its traditional application, which uses periodic 
boundary conditions. 

We end this section by noting that solving Eq. (|21|l for 
the self-energy usually requires numerical schemes. Our 
approach is to rewrite this equation in the form 



S(c^) 



dxVix) 



j <lxV{x)- 



1 - (x- Y.{uj))g(i{u)) 
1 



(24) 



{x - T.{uj))ga{uj) 



which may be solved iteratively, using as starting value 
Eo(w) = 0, the value of T,{uj) for the homogeneous sys- 
tem. The solution E(a;) to this equation in combination 
with the results described below Eq. (|20|) determines our 
CPA results for the spectra and density of states. 



IV. NUMERICAL RESULTS FOR THE 
SPECTRA 



In this paper, we restrict ourselves to the application 
of our formalism to the cylindrical aggregates occurring 
in the chlorosomes of the bacterium Chloroflexus auran- 
The application to the bilayer synthetic ag- 
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gregates of carbocyanine dyes is discussed elsewhere^ In 
terms of our stack of rings representation, the chloroso- 
mal cylinders, consisting of tens of thousands of bacte- 
riochlorophyll c molecules, have the following model pa- 
rameters (see Refs. IT^Ii^ and our discussion in Ref. 0): 
N2 = 6, h = 0.216 nm, a = 189.6°, /3 = 36.7°, and 
7 — 20°. The radius is given by i? = 2.297 nm, while the 
length may extend to hundreds of nanometers. Finally, 
we use an average single-molecule transition frequency loq 
that agrees with a wavelength of 660 nm^^i^^i^^i^ while 
for the dipole squared of a single molecule we have used 
/i^ « 20 Debye^ . The latter value was obtained from the 
integrated extinction coefficient of monomeric solutions 
of BChl c^^siLiS using the expression from Ref. lisl 

Figure 2 presents the optical spectra for several cylin- 
der lengths A^i calculated using these model parameters 
and converted to a wavelength scale to facilitate the com- 
parison to experiment 42i All intermolecular dipole-dipole 
interactions were accounted for in these spectra. The 
solid lines represent the results obtained by numerical 
simulation, where we averaged over 1000 disorder real- 
izations and used the rigorous smoothening technique 
proposed by Makhov et alm^^ to reduce the noise in 
the spectra. For the disorder strength we used a = 600 
cm~^, a value that was chosen such that the theoretical 
absorption lincwidth agrees with the one observed exper- 
imentally at room temperature. Also plotted in Fig. 2 
(dashed lines) arc the results obtained using the CPA 
modified for open boundary conditions, as described in 
Sec. mil Finally, the dotted lines represent the spectra 
obtained within the traditional CPA, which assumes pe- 
riodic boundary conditions in the ni direction (at finite 
A'^i values). 

As is clear from Fig. 2, the modified CPA is in excellent 
agreement with the numerical simulations for all spectra 
and at all sizes. Both shape and position of the spectral 
features agree almost perfectly. We conclude that the 
finite-size effects in the spectra are captured in an excel- 
lent way by the approximation made in Eq. (|23|l . These 
results justify the use of the CPA to model the spectra 
of wider cylinders^ where the size of the one-exciton 
space becomes so large that the brute-force numerical 
simulation of the spectra becomes computationally too 
expensive. 

To demonstrate that the proposed modification of the 
CPA is in fact essential to cover the finite-size effects, 
we have also applied the CPA in the traditional way, im- 
posing periodic boundary conditions in the ni direction. 
In an attempt to still account for finite-size effects, we 
have however kept the length A^i of the cylinder finite. 
We thus used Bloch waves for the longitudinal wave func- 
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Wavelength (nm) Wavelength (nm) Wavelength (nm) 

FIG. 2: Absorption, LD, and CD spectra calculated for cylindrical aggregates with the geometry of the rod elements in 
the chlorosomes of Chlorofexus aurantiacus (see text for details) and a disorder strength of cr = 600 cm~^. The solid lines 
are obtained by numerical simulations (average over 1000 disorder realizations) while the dashed lines are obtained using the 
modified CPA discussed in Sec. If If I The dotted lines represent the usual implementation of the CPA, which starts from periodic 
boundary conditions. From top to bottom, the cylinder length is varied as follows: A'^i — 15, 50, 85, 150, and 250. 



tions (fiki {ni; fc2)^The results for the absorption and LD 
spectra are shown as dotted curves in Fig. 2. For the CD 
spectrum, the use of periodic boundary conditions may 
be shown to be consistent only in the limit iVi — > oofi^ 
which is why dotted curves are absent in the CD panels. 
It is clear from Fig. 2 that the CPA with periodic bound- 
ary conditions yields spectra that differ notably from the 
exact ones, both in the position and shape of spectral 
features. The agreement is especially bad for short cylin- 
ders, where finite-size effects are most prominent. At all 
sizes considered (even for A^i = 250) the modified CPA 
constitutes a better approximation than the one with pe- 
riodic boundary conditions. Wc note that the spectra 
obtained by using periodic boundary conditions are al- 
ways red-shifted relative to the exact and the modificd- 



CPA spectra. This results from the fact that for peri- 
odic boundary conditions every molecule interacts with 
other molecules that are at most half a cylinder length 
away, while using the correct open boundary conditions, 
the molecules near one edge of the cylinder (rii small) 
have much weaker interactions with the molecules at the 
other edge {ni w Ni). Thus, using periodic boundary 
conditions, one overestimates the effect of the interac- 
tions, which shifts the spectrum too far away from the 
monomer transition. 

As final issue, we discuss the comparison between the 
exact spectra and the experimental ones. In our previous 
workii we made this comparison neglecting the effects of 
disorder in the model. One may consider those previous 
results as one extreme case, where all spectral broadening 
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is assumed to arise from homogeneous broadening, while 
the current spectra constitute the other extreme situa- 
tion, where the widths are assumed to purely result from 
inhomogeneity. Like in the homogeneous case reported in 
Ref.lH we see that in particular the CD spectra exhibit a 
strong dependence on the length A^i of the cylinder. This 
length dependence, specifically the change of the CD line 
shape around Ni = 100, was an important point in our 
previous work, as it suggests that the strong variation in 
the reported CD spectra of chlorosomesSSi^^iSi^ re- 
sults from the fact that different samples contain chloro- 
somes of different length. Apparently this conclusion sur- 
vives the incorporation of disorder and the concomitant 
localization of the exciton states. 

We stress that at the disorder value considered the typ- 
ical exciton localization size in the region of the absorp- 
tion band is several tens of molecules (Sec. 01. Thus, one 
would expect the spectra to be size-saturated at cylinder 
lengths of at most several tens of rings. For the absorp- 
tion and LD spectra, this indeed is the case, except that 
small shifts of the entire line shape still occur for longer 
cylinders as a result of the long-range dipole-dipole inter- 
actions. The slower size saturation for the CD spectrum 
results from two aspects. First, the presence of the inter- 
molecular distances in the expression for the rotational 
strengths (cf. Eq. Hll|l ') contribute to a prolonged size 
dependence. Second, being a difference spectrum, the 
CD spectrum is much more sensitive to the already men- 
tioned small shifts in the exciton energies that result from 
the long-range interactions. The fact that for clorosomes 
the length dependence of the CD spectrum survives the 
inclusion of disorder was also suggested by Prokhorenko 
et aliiS They based this conclusion on a study of the so- 
called CDM-matrix over a limited length interval, rather 
than a direct study of the spectrum. 

While the spectra for the model with disorder follow 
the same general trends as those without disorder, the 
more detailed comparison to experiment is better for the 
case with disorder. First, in the presence of disorder 
the high-energy dip in the CD spectra for A^i > 100 is 
seen to have a smaller amplitude than the low-energy dip, 
while this ratio is opposite for the homogeneous caseiii 
In experiment, the type of CD spectra with two nega- 
tive dips indeed always have a smaller amplitude for the 
high-energy dip. The effect of disorder is to smear out 
the high-energy dip, giving it a smaller amplitude. Sec- 
ond, the presence of disorder gives the absorption and 
LD spectra a more pronounced high-energy tail than is 
obtained for the homogeneous model. Indeed, these tails, 
quite typical for disordered J aggregates, are observed in 
experiment. 

To finish this section, we present in Fig. 3 the den- 
sity of states for the disordered model (solid line, ob- 
tained from simulations) and the homogeneous one (dot- 
ted line), both for A^i = 250. As above, the disorder 
was taken to be cr = 600 cm~^; the sticks in the den- 
sity of states for the homogeneous case were convoluted 
with Lorentzian curves of FWIIM= 20cm~^. The dis- 
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FIG. 3: Density of states for homogeneous (dashed line) and 
disordered (solid line) cylindrical aggregates of the chlorosome 
structure with a length of A*'i — 250 rings. In the former case, 
the excitons were homogeneously broadened by a Lorentzian 
of FWHM 20 cm^ ^ , while in the latter case a disorder strength 
of cr = 600 cm^^ was used. The disorder results were obtained 
from numerical simulations, averaging over 1000 realizations. 

order clearly smears the discrete peaks still visible for 
the homogeneous case and leads to one broad feature 
that peaks somewhere in the middle of the exciton band. 
This is in marked contrast to the frequently studied one- 
dimensional aggregates, where the density of states, even 
in the presence of disorder, peaks at the band edges. So, 
even while the cylinder constitutes a strongly anisotropic 
system, its density of states deviates from a simple one- 
dimensional picture. 



V. NUMERICAL RESULTS FOR THE 
LOCALIZATION CHARACTERISTICS 

We now turn in more detail to the localization of the 
exciton wave functions. As a first step, we consider the 
participation ratio >C~^(w). In Fig. 4, we have plot- 
ted 9C~^{uj)/A for chlorosomes within the homogeneous 
model (Fig. 4(a), A^i = 250) and the disordered model 
(Fig. 4(b), A^i = 150, 200, 250, and 300). In the homoge- 
neous model we replaced the delta functions in Eq. itT^ 
and the density of states by Lorentzians with a FWHM of 
20cm~^, while in the disordered case (disorder strength 
a — 600 cm~^, as in the previous section), we used the 
smoothening tcchnique^Si^ to reduce the noise in the 
simulations. The normalization factor 9/4 was intro- 
duced to guarantee that for the homogeneous case we 
recover the total number of molecules N in the cylinder, 
see discussion below Eq. (I12II . This is clear from Fig. 4(a), 
where, indeed, inside the exciton band 9£~^(w)/4 ob- 
tains an almost constant value of 1500 molecules, the 
total number of molecules in a cylinder of 250 rings. 

In Fig. 4(b) we see that the disorder strength of 
600 cm~^ leads to clear localization of the exciton states 
at all energies, even in the centre of the band, where the 
robustness against localization always is strongest. The 
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FIG. 4: (a) Energy dependent participation ratio for homoge- 
neous cylindrical aggregates of the chlorosome structure with 
a length of A^i — 250 rings. The dotted line indicates the total 
number of molecules (1500) in the cylinder, (b) As in panel 
(a), but now in the presence of diagonal disorder of strength 
a = 600 cm~^ and considering four cylinder lengths. From 
top to bottom the curves correspond to A^i = 300, 250, 200, 
and 150, respectively. 



fact that the band edge states are rather strongly local- 
ized is not surprising: the ratio of the disorder strength 
a and the total exciton bandwidth is approximately 0.26. 
In the spectral region where the absorption band occurs 
(720 — 750 nm), the participation ratio yields a number 
of several tens of molecules over which the exciton states 
are delocalizcd. We emphasize that this number varies 
substantially (by a factor of 8) over the width of the ab- 
sorption band. On the other hand, it is seen that in this 
energy region, the participation ratio hardly depends on 
the cylinder length anymore, which is in accordance with 
the fact that the exciton states are localized on a region 
of the cylinder that is much smaller than its total size. At 
the peak of the absorption spectrum (about 740 nm), the 
calculated participation ratio implies that the excitons 
are shared coherently by about 15 molecules (the auto- 
correlation function yields 11 molecules, see below). This 
value may be compared with the number of 7.4 — 7.6 ob- 
tained from measuring the bleaching ratio in absorption 
difference experiments Given the large variation in 
the localization size over the absorption band, the arbi- 
trariness and uncertainty present in any definition of a 
localization size, and the nonlinear nature of the absorp- 
tion difference experiments, the agreement between our 



simulations and these experimental data is good. Our 
localization size is considerably smaller than the 40 — 50 
molecules obtained in the simulations of Prokhorenko et 
ali^ Still, for the same reasons as stressed above already, 
considerable room exists for deviations in the values re- 
ported from different model studies. 

To obtain insight into the possible anisotropy of the 
localization properties, we also studied the autocorrela- 
tion function C{n;uj) defined in Eq. H13|) . We note that 
this anisotropy may also be studied by plotting individ- 
ual wave functions, as we did in Fig. 10 of Ref. Iisl but 
such plots have the drawback that they represent arbi- 
trarily picked and hopefully typical states, while the cor- 
relation function gives statistical information. In Fig. 5 
we present three-dimensional plots ((a)-(c)) as well as 
contour plots ((d)-(f)) of C(n;w) for frequencies corre- 
sponding to wavelengths of, respectively, 700 nm, 740 nm, 
and 780 nm for a cylinder of length A^i = 250 and a dis- 
order strength of cr = 600 cm^^. To make these plots, 
the cylinder was cut along a line parallel to the z axis 
(the ni axis in the plots) and unwrapped. 

Several observations can immediately be made from 
Fig. 5. First, with increasing energy the wave functions 
clearly get more extended, which is in agreement with 
Fig. 4(b). In fact, if we calculate N2,i{uj) (defined at the 
end of Section ini) from C{n;uj), we find for A^i ~ 250 
the values iV£,(700nm) = 147, A^£,(740nm) = 11, and 
iV^j(780nm) = 3. These values are in reasonably good 
agreement with those obtained from the participation ra- 
tio: 9/:-i(700nm)/4 = 119, 9/:-i(740 nm)/4 = 15, and 
9/:-^(780nm)/4 = 2.9. We also note that the values for 
Ngiiuj) (and the plots for C(n; w)) at 740 and 780 nm do 
not depend on the cylinder length anymore at A'^i = 250. 
Second, as is most clearly visible in the contour plots, 
at the short length scale the wave functions exhibit a 
clear anisotropy along a direction that is neither given 
by ni=constant, nor by n2=constant, i.e., neither in the 
ring direction, nor in the direction of the helices drawn 
as dashed line in Fig. 1. In fact, the slanting direction 
of the contours relative to the vertical axis observed in 
Fig. 5 (which varies slightly with varying energy) , closely 
resembles the slanting of the equal-phase lines of the wave 
functions in the ^2 = ±1 bands of the homogeneous cylin- 
der discovered in Ref.0. In the remainder of this section, 
we will explain this behavior, starting with a perturba- 
tive picture in which weak disorder mixes exciton states 
within and between the ^2 bands that exist for homoge- 
neous cylinders. 

In Ref. [T3. we have shown that the optically dominant 
exciton states in the fc2 = and fc2 = ±1 bands of the ho- 
mogeneous cylinder are well approximated by analytical 
expressions of the form Eq. H20|l with the ansatz 



iVi + 1 + 1 



(25) 

Here, fci is a positive integer, small compared to A'^i for 
the states of interest, s = ^2/1^2!, 7 is the helical angle of 
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FIG. 5: Three-dimensional plots ((a)-(c)) and contour plots ((d)-(f)) of the autocorrelation function C{n;uj) [Eq. CHU for 
cylindrical aggregates of the chlorosome structure with a length of A'^i = 250 rings and a disorder strength of a = 600 cm^^ 
(averaged over 150 realizations) at three different energies, corresponding to 700 nm ((a) and (d)), 740 nm ((b) and (e)), and 
780 nm ((c) and (f)). The cylinder surface is represented by cutting it along a line parallel to the z axis and unwrapping it. 
This cutting line is the ni axis in the plots; points with constant value of ni lie on the same ring of the cylinder. Due to the 
helical structure of chlorosomes, lines of constant n2 make a finite angle with the ni direction (compare the dashed line in 
Fig. 1). The contour plots distinguish 25 equally large intervals for the value of C{n;ui), represented on a grayscale (see legend 
above the contour plots). 



the cylinder (see Fig. 1), and is a phase angle that is 
used to optimize the quality of the ansatz for the complex 
effective one-dimensional Hamiltonian of the band with 
wave number The ansatz Eq. (|25() is exact if trans- 
fer interactions only occur between molecules on nearest- 
neighbor rings; in that case, the total angle + \k2\l is 
the phase of the effective nearest-ring interaction. Opti- 
mizing approximately accounts for mixing of the wave 
functions Eq. (|25|l by non-nearcst-neighbor interactions. 



The angle 9k2 was found to be responsible for the slanting 
of equal-phase lines of the wave functions in the fc2 = ±1 
band of the homogeneous cylinder (chirality of the wave 
functions) relative to the vertical axis, which seems es- 
sential to understand the behavior of C(n; uj) observed in 
Fig. 5(d)-(f). We therefore reconsider the ansatz solu- 
tions Eq. (|25|) , with a special focus on 9k2 and extending 
the treatment to the bottom regions of all k2 bands. 

The k2 = band. The Hamiltonian for the k2 — 
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FIG. 6: Gray-scale density plots of the squared amplitudes 
(darker represents a higher amplitude) of the lowest states 
in the six exciton bands of homogeneous cylindrical aggre- 
gates of the chlorosome structure with a length of A^i — 250 
rings: ^2 = (a), k2 = ±1 (b), k2 = ±2 (c), and k2 — S 
(d). Plots are given on the unwrapped cylinder surface. The 
wave functions were obtained from exact diagonalization of 
the corresponding homogeneous Hamiltonians and are well 
described by Eq. I2UII . with the ansatz Eq. I25II (see text for 
details). Relative to the average molecular transition energy 
ujo, the lowest states have the energies — 1288.73cm~^ (a), 



-1291.92 cm-i (b), 
(d). 



-1244.97 cm"^ (c), and -1166.09 cm- 



band is real, implying that also the wave functions can 
be chosen real. Hence, we have 9q ~ 0. It was shown in 
Ref. that the bottom state of this band (the fci = 1 
state, which carries 81% of the oscillator strength in this 
band) is well described by the ansatz. For this band, the 
wave function has equal amplitude on all molecules of 
a certain ring, while the quantum number ki gives the 
number of maxima along the z axis of the cylinder. This 
is illustrated in Fig. 6(a), where we plotted the amplitude 
squared of the exact wave function (obtained by numer- 
ical diagonalization) of the lowest state in the fc2 = 
band of the homogeneous cylinder. 

The A:2 = il bands. The effective Hamiltonian for 
these two degenerate bands is essentially complex and 
the phase factor in Eq. (|25() is needed. In Ref. we 
have shown that the superradiant states (three in each 
band) all occur near ki = k* = nint[|6'i|(A^i -I- l)/7r] (nint 
denoting the nearest integer function), which is a finite 
energy above the bottom of the ki = ±1 bands. Excellent 
agreement was found between the ansatz wave functions 
and those for the superradiant states of the homogeneous 
cylinder obtained by numerical diagonalization for the 



optimized value of 6±i = 4.3° Presently, we are not 
only interested in the superradiant states, but even more 
so in the states near the bottom of the ^2 = ±1 bands 
(fci in the order of unity), because these states are closest 
to the bottom of the ^2 = band and should be expected 
to mix with those states in the presence of (weak) dis- 
order. We may apply the same approach as followed in 
Ref. for the superradiant states to the states near the 
band bottom. Thus, we fix 9±i such that the mixing be- 
tween the bottom states in the fc2 = ±1 bands resulting 
from the long-range interactions does not diverge. This 
eventually yields 6'^^ = 3.8° (for Ni = 250), which only 
slightly differs from the phase angle in the superradiant 
region. 

If we plot the amplitude squared of the exciton wave 
function on the unwrapped cylinder, a state of the form 
Eq. H2U|) with Eq. H25|) gives lines of equal intensity that 
arc dictated by lines of equal phase 3(0^3 -I- |fc2|7)ni -I- 
27rA:2"-2/^2- The angle x of these lines with the z axis is 
easily derived to obey 



tanx 



(26) 



This implies that the wave function rotates around the 
cylinder over a number rings given by 



(27) 



From this we find that in fc2 = ±1 bands w 84 for the 
wave functions of the superradiant states and n{ w 95 for 
the bottom states. In Fig. 6(b) we show the amplitude 
squared for the bottom state of the k2 = I band obtained 
by numerical diagonalization. Of course, on every ring we 
find a modulation with two maxima, as is appropriate for 
this band. More importantly, we find good agreement 
between the value of that may be obtained from this 
plot and the above estimate obtained from the ansatz 
wave function. 

The k2 = ±2 bands. These bands, which contain no 
oscillator strength in the homogeneous limit, were not 
considered in Ref. El They may be treated in a way 
similar to the fc2 = il bands. If we focus on optimizing 
the ansatz at the bottom of the bands, we find 6±2 = 
6.6° for A^i = 250 and that for this value the ansatz 
Eq. for fci = 1 indeed gives a good description of 
the numerically obtained lowest state. From Eq. H27|) 
and 9±2 = 6.6° a value of nl w 109 can be estimated, 
which is seen to be in good agreement with the exact 
wave function plotted in Fig. 6(c). Of course, the states 
in the k2 = ±2 bands exhibit a modulation of squared 
amplitudes inside each ring with four maxima. 

The ^2 = 3 band. This band, too, has no oscillator 
strength. Like the fc2 = band, the fc2 = 3 band is non- 
degenerate and is governed by a real Hamiltonian. Thus, 
^3 + 37 = and, like for the fc2 = case, no optimiza- 
tion condition for ^3 has to be solved. Figure 6(d) gives 
the amplitude squared of the numerically obtained lowest 
state in this band, which exhibits no modulation inside 
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FIG. 7: As in Fig. 6, but now including diagonal disorder 
of a strength a — 20 cm^^. The weakness of the disorder 
allows for an identification to the k states of the homogeneous 
aggregate in Fig. 6. The energies of the plotted states are 
-1288.96 cm"^ (a), -1292.67 cm"^ (b), -1245.97 cm"^ (c), 
and -1167.55 cm-i (d). 



rings (the Bloch factor is (—1)"^ which gives unity upon 
taking the square). One does observe a vertical modula- 
tion with a high periodicity, which is due to the fact that 
for this band the lowest state is not the ki ~ 1 state, 
but rather a state that corresponds to a high value of fci . 
This is due to the fact that for this band the effective 
Hamiltonian has dominant positive in stead of negative 
interactions. 

The important conclusion of the above analysis is that 
the bottom states in the k2 = and fc2 = 3 bands ex- 
hibit no chiral behavior (no equal-intensity lines slanted 
relative to the z axis), while the bottom states of the 
A:2 = ±1 and the ^2 = ±2 bands do show a chirality 
which, moreover, is very similar in magnitude as esti- 
mated by the number n^. If we now allow for weak dis- 
order, mixing of the k states within and between bands 
will occur. Focusing on the bottom region of the den- 
sity of states, we thus find that the wave functions in the 
homogeneous ^2 = and fc2 = 3 bands will (perturba- 
tively) acquire a slanted contribution as well. In view of 
the fact that the slant angle varies little inside as well as 
between the A:2 = ±1 and k2 = ±2 bands, we thus ex- 
pect a similar slanting for all low-energy states at weak 
disorder. This is clearly confirmed by Fig. 7, where we 
presented the bottom states of the k2 = 0, 1, 2, and 3 
bands weakly perturbed by the presence of a very small 
disorder strength of cr = 20 cm~^. If we further increase 
the disorder strength, the mixing will become nonpertur- 



bative and will also involve higher-lying states, described 
by different x values and a different number of nodes 
in the vertical direction. These are the ingredients for 
the localization visible in Fig. 5; however, the slanting of 
the autocorrelation function still clearly reflects the fact 
that locally the wave functions have a chirality that de- 
rives from the homogeneous states. The reason is that 
this chirality is mainly driven by the dominant transfer 
interactions, which occur over a few rings. 



VI. CONCLUSIONS 

In this paper we have theoretically investigated the ef- 
fects of disorder on the linear optical properties and the 
localization behavior of the exciton states of cylindrical 
molecular aggregates. As specific example, we have used 
the structure for the cylindrical aggregates found in the 
light harvesting systems (chlorosomes) of green bacte- 
ria. We have calculated the absorption, LD, and CD 
spectra for various cylinder lengths in the presence of 
Gaussian diagonal disorder, using both numerical simu- 
lations and the CPA. To this end, we modified the usual 
implementation of the CPA to account for finite-size ef- 
fects (open boundary conditions). By comparison to the 
simulation results, we have shown that this new imple- 
mentation yields an excellent approximation, in fact sig- 
nificantly better than the usual CPA, which starts from 
periodic boundary conditions. We have shown that for 
the chlorosomes the inclusion of disorder improves the 
comparison to experiment in several respects. Most im- 
portantly, however, the inclusion of disorder does not af- 
fect the main conclusion of Ref.^^that the CD spectrum 
has a strong size dependence up to cylinders of hundreds 
of rings long. We demonstrated and discussed that this 
effect occurs in spite of the fact that the exciton localiza- 
tion size is far smaller than the cylinder size. 

The localization behavior of the excitons was studied 
by examining two quantities: the participation ratio and 
an autocorrelation function of the exciton wave functions. 
While the former only yields an approximate measure 
for the number of molecules participating in the excita- 
tion, the latter gives additional information about the 
direction of localization on the cylinder. In the case of 
chlorosomes, we found that the excitons that dominate 
the optical properties (the ones near the band bottom) 
have a strongly anisotropic localization behavior, being 
extended mainly along helices, whose direction is dictated 
by the interplay of the various intermolecular excitation 
transfer interactions in the cylinder. As we have demon- 
strated, this chiral behavior finds its roots in the chirality 
of the fully extended exciton states for the homogeneous 
cylinder. Both the participation ratio and the autocor- 
relation function show that, starting from the bottom of 
the exciton band, the exciton states become more ex- 
tended with increasing energy (both quantities give com- 
parable values for the energy dependent localization size). 
Also the chirality of the wave functions is energy depen- 
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dent, but this effect is small. It is of interest to speculate 
whether the chiral nature of the wave functions could be 
detected in terms of a rotating polarization of the light 
emitted by an exciton while propagating. This would re- 
quire single-aggregate experiments with high spatial res- 
olution and a time-resolved detectioni^2iSli 
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